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ABSTRACT 

We use results from a constrained, cosmological MHD simulation of the Local Universe to 
predict radio haloes and their evolution for a volume limited set of galaxy clusters and com- 
pare to current observations. The simulated magnetic field inside the clusters is a result of tur- 
bulent amplification within them, with the magnetic seed originating from star-burst driven, 
galactic outflows. We evaluate three models, where we choose different normalisations for the 
cosmic ray proton population within clusters. Similar to our previous analysis of the Coma 



cluster (Donnert et al.|2010), the radial profile and the morphological properties of observed 



radio halos can not be reproduced, even with a radially increasing energy fraction within the 
cosmic ray proton population. Scaling relations between X-ray luminosity and radio power 
can be reproduced by all models, however all models fail in the prediction of clusters with no 
radio emission. Also the evolutionary tracks of our largest clusters in all models fail to repro- 
duce the observed bi-modality in radio luminosity. This provides additional evidence that the 
framework of hadronic, secondary models is disfavored to reproduce the large scale diffuse 
radio emission, of galaxy clusters. We also provide predictions for the unavoidable emission 
of 7-rays from the hadronic models for the full cluster set. None of such secondary models 
is yet excluded by the observed limits in 7-ray emission, emphasizing that large scale diffuse 
radio emission is a powerful tool to constrain the amount of cosmic ray protons in galaxy 
clusters. 
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1 INTRODUCTION 

The thermal gas, that is the dominant component in the Inter- 
Galactic-Medium (IGM), is mixed with magnetic fields and rela- 
tivistic particles, as proven by radio observations which detected 
Mpc-sized diffuse radio emission from the IGM, in the form of ra- 
dio halos and relics (e.g. Feretti 2003, Fer rari et al.|20 08 ). These 
Mpc-scale radio sources are found in a fraction of massive clus- 
ters with complex dynamics, which suggests a connection between 
non-thermal emission and cluster mergers (e.g.JBuote 200l] [Verr| 
|turi et al.|2008||Brunetti et al.|2009| > Cluster mergers are the most 
energetic events in the universe and a fraction of the energy dis- 
sipated during these mergers may be channelled into the amplifi- 
cation of the magnetic fields (e.g. |Dolag et al. 2002 Subra manian| 
|et al.|2 006; Ryu et al. 20081 and into the acceleration of relativistic, 
primary, electrons and protons via shocks and turbulence (e.g. En- 



sslin et al. 1998; Sarazin 1999; Brunetti et al. 2001 , 2004 , Petrosian 
200T]|Gabici & Blasi|20"03| |Ryu et al.|2003||Cassano & Brunetti] 



2005 ; Pfrommer et al. 2006; Brunetti & Lazarian 2007 ; Vazza et al. 



2009l > 



Relativistic protons in the IGM have long life-times and re- 
main confined within galaxy clusters for a Hubble time (e.g. Volk 
[eraL][T996l |Berezinsky etlr^[T997) ). As a consequence they are 
expected to be the dominant non-thermal particle component. Col- 
lisions between these relativistic protons and the thermal protons in 
the IGM generate secondary particles that combined with the pri- 
mary relativistic particles are expected to produce a complex emis- 
sion spectrum from radio to 7-rays (e.g. Blasi 2001 ; Brunetti et al. 
|2009[ >. Only upper limits to the 7-ray emission from galaxy clus- 
ters have been obtained so far (Reimer et al. 2003; Perkins & the 
IVERITAS Collaboration||2006l [Aharonian 2009bja] |The MAGIC" 
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Collaboration & Aleksic 2009) however the FERMI Gamma-ray 
telescope will shortly allow a step forward, having a chance to 
obtain first detections of galaxy clusters or to put stringent con- 
straints on the energy density of the relativistic protons. Most im- 
portantly, in a few years the Low Frequency Array (LOFAR) and 
the Long Wavelength Array (LWA) will observe galaxy clusters at 
low radio frequencies with the potential to discover the bulk of the 
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cluster-scale synchrotron emission in the Universe (e.g. EnBlin & 
|R6ttgering|2002l|Cassano et al.|2009||2006l >. 

The emerging theoretical picture is very complex and modern 
numerical simulations provide an efficient way to obtain detailed 
models of non thermal emission from galaxy clusters to compare 
with present and future observations. Advances in this respect have 
been recently obtained by including aspects of cosmic-ray physics 
into cosmological Lagrangian simulations mostly focussing on the 
acceleration of relativistic particles at shocks and on the relative 
production of secondary electrons (e.g. Pfrommer et al. 2008}. In 
this work we investigate the non-thermal emission from secondary 
particles in galaxy clusters extracted from Lagrangian cosmolog- 
ical simulations and, for the first time, we report on an adequate 
comparison between our expectations and observations. 



of the resolution, to cover for long range gravitational tidal forces 
arising from the low-frequency constrains. 

In the evolved density field, many locally observed galaxy clusters 
can be identified by position and mass. Especially the Coma cluster 
(see |Donnert et al.|2010) shows remarkable similarities in morphol- 
ogy. A fly-through of the simulation can be downloaded from the 
MPA WebsittQ 

The initial conditions were extended to include gas by splitting dark 
matter particles in the high resolution region into gas and dark mat- 
ter particles of masses 0.69 x 10 9 M and 4.4 x 10 9 M respec- 
tively. Therefore the biggest clusters are resolved by about a million 
particles. The gravitational softening length was set to 10 kpc. This 
is comparable to the inter-particle separation found in the centre of 
the largest clusters. 



2 SIMULATIONS 



2.2 Magnetic Fields from Galactic Outflows 



The simulation was done using the cosmological simulation code 
GADGET-2 ( |Springel|2 005 ) with a treatment for magnetic fields. It 
features an entropy conserving formulation of Smooth Particle Hy- 
drodynamics (SPH) (Springel & Hernquist 20021, which is supple- 
mented with the formulation of ideal MHD presented in |Dolag"&| 
Stasyszyn ( 2008 ). The implementation follows the induction equa- 
tion and computes the back reaction of the magnetic field using a 
symmetric formulation of the Lorentz force. We used a divergence 
cleaning scheme presented in |B0rve et al.| l |2~001) , which reduces 
numerical noise in shocks by subtracting the magnetic force which 
is proportional to the divergence of the field. It also helps to sup- 
press the clumping instability particle based MHD codes encounter 
in regions with small plasma f3 (i.e. where magnetic pressure con- 
siderably exceeds thermal pressure). 

In non radiative simulations like ours, regions with small 
plasma /3 are rare. Only the strong shocks in cores of galaxy clus- 
ters during major mergers produce enough compression to amplify 
the field to become dynamically dominant. These mergers are rel- 
atively brief events and are handled more accurately with our new 
numerical treatment (see Dolag & Stasyszyn 2008 , for details). 
Borga nT~et al.| {2006[ > have shown that non radiative simulations 
overpredict the gas density in cores of galaxy clusters. This af- 
fects our simulation as well and can be seen in density and mag- 
netic field profiles as well as in X-ray luminosities. As cosmologi- 
cal MHD SPH simulations lack physical dissipation, radiative SPH 
MHD simulations are not feasible at the moment. On the other hand 
secondary models have difficulties reproducing the outer parts of 
radio haloes correctly. Therefore our main focus lies on these re- 
gions, where the simulations are not affected by the overpredicted 
gas density. 



2.1 Initial Conditions 



We used a constrained realisation of the local universe (see |Dolag| 
|et al. | j2005 1 and references therein). The initial conditions are sim- 
ilar to those used in Mathis et al. ( 2002) to study the formation of 
the local galaxy population. They were obtained based on the IRAS 
1.2-Jy galaxy survey. Its density field was smoothed on a scale of 
7 Mpc, evolved back in time to z — 50 using the Zeldovich approx- 
imation and assumed to be Gaussian (Hoffman & Ribak 1991}. The 
IRAS observations constrain a volume of w 115 Mpc centered on 
the Milky Way. It was sampled with dark matter particles and em- 
bedded in a periodic box of as 343 Mpc comoving. Outside of the 
inner region, the box is filled with dark matter particles with 1 /6th 



The origins of magnetic fields in galaxy clusters are still under de- 
bate. It is assumed that some kind of early seed magnetic field is 
amplified by structure formation through adiabatic compression, 
turbulence and shear flows to values observed today (~ 1 — 10 /iG 
in clusters). Three main classes of models for the seed field exist: 
At first the seed fields can be created in shocks through the "Bier- 
mann battery " l |Kulsrud et al.|1997||Ryu et al.|19"98||Miniati et al.| 
|2001| >. A second class of models invokes primordial processes to 
predict seed fields that fill the entire volume of the universe. The 
coherence length of these fields strongly depends on the details of 
the model (see |Grasso & Rubinstein|200l| for a rev iew). Finally 



the seed can be produced by AGN ( EnBlin et al. 1997 Furlanetto & 
|Loeb|200"T) or starbursting galaxies ( TVolk & Atoyan 2000) at high 
redshift (z w 4 — 6), whose outflows contaminate the proto-cluster 
region. 

Cosmological simulations using SPH ( Dolag et aL|1999[|2002] 
2005) and grid based Adaptive Mesh Refinement (AMR) codes 



(Brug gen et al.|2005[|Dubois & Teyssier|2008||Li et al.|2008| > were 
able to show that observed Faraday rotations are compatible with 
a cosmological seed field of « 10~ n G. They also suggest that 
spatial distribution and structure of cluster magnetic fields are de- 
termined by the dynamics in the velocity field caused by structure 
formation (Dolag et al.|1999| 



For this work we follow 



2002) 



Jonnert et al. ( 2008 1 in terms of mag- 



netic field origin. They use a semianalytic model for galactic winds 
( |Bertone eT al. 2006) to seed magnetic fields in a constrained cos- 
mological MHD SPH simulation. The continuous seeding process 
is approximated with an instantaneous seed at z ~ 4. As they were 
able to show, the main properties of magnetic fields obtained in 
clusters were not influenced by that approximation. 
The wind model used assumes adiabatic expansion of a spherical 
gas bubble with homogeneous magnetic energy density around ev- 
ery galaxy below a certain mass threshold. The magnetic bubble 
can be characterised by radius and field strength. The galaxy in- 
jects gas into the bubble earring frozen-in magnetic field from the 
disc into the bubble over the star-burst timescale. Its final size is 
determined by the wind velocity, which is a function of the star 
formation rate and the properties of the ISM. |Bertone et alT| ( 2006) 
give an evolution equation for the magnetic energy in the bubble de- 
pending on the star-burst timescale. The energy is converted into a 
dipole moment and seeded once at a chosen redshift. The magnetic 



1 http://www.mpa-garching.mpg. de/galform/data_vis/index.shtml#moviel2 
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Figure 1. Full sky maps of the simulation in galactic cordinates. From top to bottom, electron density, temperature and magnetic field, projected through the 
whole box. The inlay shows a zoom onto a 3x3 degree region arround the Coma cluster, respectively. In the upper most map, the most prominent clusters of 
the local universe are labled and the arrow in the inlay points towards north. 
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Parameter 


Value 


Source 


Ro 


400 pc 


(Klein et al.|l988j 


Bo 


5mG 


{Donnert etal.|2008 




B G 


3mG 


Jponnert et al. 2008 




M* 


10M Q /yr 


jde Griis|200li 




*sb 


150 Myr 


jde Grijs 2001 


MlSM 


< 300 X 1O 12 M 


from simulation 



Table 1. Summary of the parameters used for the wind model. This corre- 
sponds to the 0.1 Dipole simulation in I Donnert et al. 2008), which fits best 
to observations of Faraday rotation 
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Figure 2. Radial profile of the density in our simulated Coma cluster. In red 
we plot the best fit beta model (po = 6.4 X 10 15 , = 0.8, r CO rc = 0.06), 
excluding the density bump at 0.6 i? v ir . A fit to observations of Coma by 
|Mohr et al.| (1999) is shown in green. Due to the non-radiativity of our 
simulation the gas density is overestimated by a factor of 5 in the center 



field is then amplified by structure formation to /j,G level. For de- 
tails on the wind model refer to Bertone et al. ( 2006 ); Donner Fet al.| 
(2008). Figure [T] shows full sky maps produced from the simula- 
tion, projecting the electron density, temperature and the magnetic 
field. The magnetic field closely follows the density distribution. 
The magnetic field is more patchy in the filaments compared to a 
cosmological seed because the seeding by individual galaxies does 
not overlap. 

The simulation used in this work is based on the 0.1 Dipole 
parameter set from Donnert et al. ( 2008), shown in table [T] It rep- 
resents the best fit to observations of Faraday rotation presented 
in |Donnert et al.|j2008[ >. There the parameter space was explored 
based on observations from the wind in M82. It was demonstrated 
that the resulting magnetic field does not critically depend on the 
exact choice of the parameters. 



MODELLING HADRONIC SECONDARY ELECTRONS 
IN CLUSTERS 



In secondary models proposed for the origin of radio halos ( |Den-| 
mson [19801 |BlasT& Colafrancesco|1999| >, the emitting cosmic ray 
electrons (CRe) are a product of cosmic ray proton collisions (CRp) 
with thermal protons in the cluster. Possible mechanisms for the in- 
jection of CRp in the IGM include shock waves caused by cluster 
accretion and merger s (e.g.|Gabici & Blasi|2003| |Ryu et al.|2003| 
|Pfrommer et al.|2 006 , Vazza et al. 2009) outflows from radio galax- 



ies (AGN) (e. g.|Ensslin et al.|1998||Rachen|2008l > and supernova 
driven winds ( |Volk et al.|1996) . The resulting CRp spectrum from 
these theoretical scenarios is a power law (e.g. Schlickeiser 2002 ), 
that is also consistent with observations of galactic cosmic rays. 

We assume a power law spectrum of CRp in the IGM and 
calculate the spectrum of high-energy secondary electrons under 
stationary conditions by considering only synchrotron and inverse 
Compton losses and by neglecting re-acceleration processes in the 
IGM (Section|3TT). 

Throughout the first two of our models, we use the magnetic fields 
obtained from the constrained cosmological simulation described 
in section \22\ In a third model we increased the magnetic field, 
especially at larger distances to the center, by up-scaling the 
magnetic field in our simulation to be B oc <J~p, which allows for 
better comparison with results obtained in previous literature, as 
there such a relation is often used (as for example in Pfrommer 
[etal1 ( |2008) ). 



3.1 Synchrotron emission from secondary models 

We consider a spectral distribution of the population of relativistic 
CRp described by a power law in energy E p : 



N(E P ) = 



K P E P ' 



(1) 



if not stated otherwise, we use a spectral index a p = 2.6. The 
normalisation K p is model dependent and is defined in the next 
Sections. 

The main channel of hadronic interaction between the CR protons 
and the ambient medium is multi-pion production jBlasi & Co-| 
lafrancesco 1999): 



Pcr + pth 
± 



, 4 
7T + 7T 



7r + anything 



e + v M + v e 
27 



In this case the injection spectral rate of secondary e is (e.j 
|Brunetti & Blasi|2005) : 



Qe(E) 



) 7V (E 7T )dE^x 



dE ll F e (E n , Efj,, E e )F IJ/ (E li , E K ), 



(2) 



where F^(E e , E^, E n ) is the spectrum of electrons and positrons 
from the decay of a muon of energy E^ produced in the decay of 
a pion with energy E^ (taken from Blasi & Colafrancesco ( 1999 ), 
Fy, (Ep , En ) is the muon spectrum generated by the decay of a pion 
of energy E^ (e.g. Moskalenko & Strong 1998 ), and the pion in- 
jection rate due to p-p collisions is (e.g. | Dermer| 19 86 ; Blasi & Co-| 
|lafrancesco|1999||Brunetti & Blasi|2005> 



n(Erc) = n th c / dE p N(E p )p p F TC (E TC ,E p ) x 



(3) 
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where n t h is the number density of thermal protons, a pp is the p- 
p cross-section, and F v is the spectrum of pions from the collision 



between a CRp of energy E p and thermal protons (e.g. Moskalenko 
& Strong|[T9^8l |Blasi & Colafrancesco|[T999"l |Brunetti & Blasi 



2005| >. 



The observed synchrotron emission in radio halos at frequencies 
of ~ GHz requires energies of the secondary particles of several 
GeV, considering field strengths of ~ /iG typical of the IGM (e.g. 
|Carilli & Taylor|2002) . In this case a scaling model is appropriate 
to describe the pion-spectrum from p-p collisions, and we calculate 



Q e (E) from Eqs.[2p1following |Brunetti & Blasi||2005| > and using 
a pion-spectrum : 



1 

2-Ett 



ci(l - ^) 3 ' 5 + c 2 exp(-18 
E, 



Ett , 

Ei, 



(4) 



where ci = 1.22 and C2=0.92 (Berezinskii & Kudriavtsev 19901. 
The resulting electron source function can be approximated by : 

Q c ± (E) = cn th K p E- a *>f{E, a p ) (5) 

where f(E,a p ) accounts for the log-scaling of the p-p cross- 
section at high energies and causes the spectral shape to be slightly 
flatter than E~ a ' p ; an analytical expression for the asymptotic form 
of / (for E^ » m^c 2 , E n » m^c 2 , E p » m p c 2 ) is derived 
in |Brunetti & Blasi|p005) . 

The steady state spectrum of high energy electrons in the IGM 
is given by (e.g. |Dolag & Enfflin|2000) : 



N C (E) = E(E)\ 



dE' Q C ±(E') 



(6) 



where the relevant cooling processes involve synchrotron and in- 
verse Compton losses: 



E(E) 



Sm 2 c 4 



£J i -"CMB 

8tt 8tt 



E z 



(7) 



B is the local magnetic field strength and Bq MB /8ty gives the en- 
ergy density of the CMB expressed as an equivalent magnetic field. 
We calculate the electron spectrum from Eqs.|2]|4]and|6}|7j this can 
be approximated by : 



Here Ks denotes the modified Bessel function of order 5/3. 

3 

For a power law spectrum of CRp, N P (E P ) oc E p p , the re- 
sulting synchrotron emission from secondary electrons is j u oc 

,-(a p -A)/2 



3.2 7-rays from hadronic interactions 

As already mentioned before, CRp - proton collisions produce neu- 
tral pions, which in turn decay to 2 photons. 7-rays are a direct 
measure of the CRp and provide a complementary constraint to the 
injection process of secondary electrons. 

In order to allow for a prompt comparison with recent results we 
follow the formalism described in Pfrommer & EnBlin (20041 to 
estimate the 7-ray flux from CR in our simulations. 
The 7-ray source function is : 



GeV 



(TppCnthKp 



3a 7 



(£p. 



E„ 



GeV 



2E-, 



+ 



2E„ 



(12) 



where a 7 ~ a p is the asymptotic slope of the 7-ray spectrum, 
which resembles the slope of the proton spectrum ( |Dermer|1986| >. 
The shape parameter, which describes the semianalytic model near 
the pion threshold, is <5 7 = 0.14a 7 16 + 0.44by using an effective 
cross-section opp = 32 x (0.96 + exp(4.4 — 2.4a 7 )) mbarn. The 
integrated 7-ray source density A 7 is then obtained by integrating 
the source function over energy ( Pfrommer & EnBlin 2004) : 



E-2 



A 7 = j di? 7 g 7 (i? 7 



CTppm^C 3 UthKp (-Ep,min) " p -Ep.min ( m n C 



3q 7 5 7 a p - 1 2 a i- 1 



GeV V GeV 



a-y + 1 a 7 
2(5- 



2<5 7 



(13) 



where B x (a,b) denotes the incomplete beta-function and 

[/(*)]? = f(a) - f(b). 



N e (E) = 



<7t 



l -n fh K 9(E,a P ) 



CMB 



where g is related to / in Eq.[5]via 
1 



g(E,a p ) 



dXX~ a vf(X,a p ) 



oc E 

and A w 0.2 for E w few GeV. 
The radio emissivity is : 

Emax I 



3» 



VSe 3 B 



II 

-Bmin 

J2, 



&E e &6 sin6F ( — ) N e {E) 



(8) 



(9) 



GO) 



where v c = (3/4ir)p 2 eB sin6 l /(mc) 3 is the critical frequency and 
F the integral over the synchrotron kernel : 



F(x) =x / Kn(£)d£. 



(11) 



3.3 The three models 

To investigate the dependence of the predicted properties of non 
thermal emission of clusters on the underlying assumptions, we in- 
vestigate 3 models for the distribution of magnetic fields and cos- 
mic rays in clusters. They are chosen to encompass the reasonable 
range suggested by theoretical and observational findings. 
We keep the spectral index fixed to a p — 2.6 in order to be able 
to match the typical spectrum of giant radio halos, a ~ 1.2 — 1.3 
(e.g. Ferrari et al. 20081, although a fraction of presently known 
halos has a steeper spectrum (e.g. |Brunetti et aLj 2008, Giov annim] 
|et al.|20 09 ) Also, a spectral index a p = 2.6 allows to fit the spectral 
shape of the Coma halo at v ^ 1.4 GHz, although also in this case 
the spectrum steepens at higher frequencies (e.g. |Thierbach et al.| 
[20031 iDonnert et al.|2010] > 



3.3.1 Model 1: Constant Xcr 

In our first model the energy density of the CR protons is taken as a 
constant fraction of the thermal energy density. This is reasonable 



6 /. Donnert, K. Dolag, R.Cassano, G. Brunetti 




-1M -Jill -19,0 -17.1 -IS.!) -13JI -1L0 



Figure 3. Full sky maps of the simulation in galactic cordinates. From top to bottom, X-ray, radio and 7-ray surface brightness are shown. For the radio and 
7-ray emission the model with the constant cosmic ray energy fraction was chosen. The 7-ray emission is evaluated for VERITAS (e.g. E > 100 GeV). The 
inlay shows a zoom onto a 3 X 3 degree region around the Coma cluster. In the upper most map, the most prominent clusters of the local universe are labled 
and the arrow in the inlay points towards north. 
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Figure 4. Energy in cosmic ray protons relative to thermal energy over clus- 
ter radius, inferred from Pfrommer et al. (2007). The dependence was de- 
rived assuming a linear relation between pressure and energy density. 



Figure 5. Scaling functions / sca i (r) of all clusters. It modifies the magnetic 
field globally, so that its magnetic energy scales linear with thermal energy. 



if a constant fraction of the energy that is channeled into the IGM 
to heat the gas goes into the acceleration of CRp. 

Therefore, in this model, the normalisation K p is chosen to 
have a constant fraction, X p = const, of kinetic CRp energy den- 
sity e p to thermal energy density e t h of the IGM: 



(£p. 



(14) 
(15) 



For the magnetic field distribution within the IGM we take 
directly the magnetic field extracted from the simulations. 



3.3.2 Model 2: Varying X C r (r) 

In a second model we adopt a radius dependent cosmic ray energy 
density fraction, -X"cr(V), as obtained from simulations of CR ac- 
celeration in structure formation shocks by Pfro mmer et al.|p007| >. 
We also assume a constant CRp spectral index over the whole clus- 
ter volume, a p = 2.6. 



|Pfrommer et al.| ((2008 l simulated the injection of CR protons 
by merger shocks during structure formation. They find that cosmic 
ray pressure increases relative to thermal pressure with increasing 
distance to the cluster centre. Assuming an ideal gas this directly 
translates into a radially increasing cosmic ray energy density frac- 
tion Xcn,(r) — eth/tp- We use these results to infer a radially 
varying normalisation, so that Eq.|14|becomes: 



X p (r)e th (r) 



K P (r) 



(-^p,min) 



(16) 
(17) 



Figure [4] shows the radial dependence of the fraction used in 
agreement with the average over two cool core clusters, as shown 
in |Pfrommer et a l. (20081. The fraction increases by more than 2 
orders of magnitude in the cluster external regions. This is expected 
to increase the extension of the diffuse radio emission generated 
by secondary electrons with respect to model 1. For the magnetic 
field distribution within the IGM we again use the magnetic field 
extracted from the simulations. 



3.3.3 Model 3: Scaling B oc ^fp 

In our third model we modify the magnetic field originally obtained 
from the cosmological simulations. 

The main reason for this is that the magnetic energy density, 
£b,o, of the IGM is often modelled to follow thermal energy den- 
sity, eth, which implies B oc ^fp. In contrast to that we find B oc p 
in the outer regions of our simulated clusters. To allow a compar- 
ison with models that assume B oc ^fp, we construct in our third 
model a function / sca i (r) which scales the simulated magnetic field 
to follow B oc yfp globally. This allows us to alter the dependence 
of the field with thermal density without loosing its inherent struc- 
ture which is a result of our simulations. 

To construct the scaling function, we fit beta models np = 
no(l + r 2 /r^ orc ) 313 ^ 2 to spherically averaged density and magnetic 
field profiles. The scaling function is then: 



/scalO) = 



(po(l + r 2 /r-c 0rc , p ) 



3/2/3 f 



1/2 



Bo(l + ryr 2 r/^ 



(18) 



and is shown in figure [5] for every cluster in the sample. This 
changes the field strength at r > 0.2R vir without strongly 
affecting the field in the innermost regions. 

For the cosmic ray distribution in this model, X p (r), we use the 
same radially dependent profile as in our model 2. 



4 APPLICATION TO THE CLUSTER SAMPLE 

Having the diffuse synchrotron emission from secondary electrons 
in our simulated clusters under control, the aim of this Section is to 
compare the simulated properties with the most relevant observa- 
tional properties of radio halos. 

We use an unprecedented sample of fairly-massive simulated clus- 
ters made by the 16 most massive objects extracted from our sim- 
ulations. This allows for a statistical comparison with the main ob- 
servational properties of radio halos obtained from studies of vol- 
ume limited samples of radio halos in galaxy clusters. Addition- 
ally our sample (as well as the observational samples) contains 
the Coma cluster, so that we can calibrate K p in all our models 
to match the observed radio luminosity of the Coma halo where 
needed. 
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Figure 6. Faraday rotation over radius for a mass selected subsample of 
our simulated clusters. Shown are the original (green) and upscaled (red) 
median over the whole sample in every radial bin. We also plot observations 
from Coma (turquoise), A119 (blue) and a sample of Abell clusters (black) 
(Kim et al.|199T]|Clarke et al.|2001||Feretti et al.|1995|| 1999| . The errors in 
Abell cluster sample were estimated using a bootstrapping technique (see 
text). 



4.1 The magnetic field in our cluster set 

The magnetic field obtained for our simulated Coma cluster is con- 
sistent with the one infered from modeling Faraday Rotation Mea- 
sures of the Coma cluster, as shown in l |Donnert et al.|2 010). 
To compare the "statistical" properties of the magnetic field in the 
complete set of simulated clusters with observations and to show 
the effect of re-scaling the field (model 3), we plot in figure [6] ra- 
dial profiles of Faraday Rotation. Here we bin all clusters with gas 
mass M > 3 x 10 14 Mq in radius and take the median in each bin. 
In green we plot the field profile from simulations, and in red the 
scaled field profile (model 3). To compare we overplot observations 
of a sample of Abell clusters from |Kim et al.H1991|l;|Clarke et all] 
^2001) , observations of A119 from |Feretti et al.| ( |1999^ and obser^ 
vations of the Coma cluster (Feretti et al. 1995). Because of the 
small number of points, the error from the Abell clusters was esti- 
mated using a bootstrapping technique. For every bin we generate 
several samples, by selecting every data point a random number of 
times and computing the median of this subsample. We then com- 
pute the standard deviation of all samples to get an estimate of the 
error in each bin. 



4.2 Radial profile of the radio emission 

Since both the target thermal-protons and the expected magnetic 
field strength in the IGM decrease with distance from the cluster 
centre, most of the synchrotron luminosity emitted by secondary 
electrons should be produced in the cluster-core region. This causes 
the radial profile of secondary-generated radio halos to be substan- 
tially steeper than those of the observed halos (e.g. |Brunetti et al.| 
|2004j >, although formally this discrepancy may be alleviated by as- 
suming that both the magnetic field and CRp have flat spatial dis- 
tributions ( |Pfrommer & Enfilin|2004) . 

|Cassano eFa l. ( 2007 ) report radial profiles of the radio emission of 
5 well studied radio halos. We convolve our synthetic radio maps 
with a Gaussian using the typical beam size (43 kpc) of the ob- 
servations, and in figure|7]compare these profiles with the observed 
ones. Both simulated and observed profiles were normalised to one. 



We also include the radial profile of the Coma halo from Deiss et al. 
(1997). 

In all cases (models 1-3), Figure|7]shows that the radial profiles of 
the secondary-generated radio halos in our simulated clusters are 
considerably steeper than the observed ones. 

In agreement with |Brunetti et aL] ( |2004) , by assuming a con- 
stant Xp (first panel), the simulated radio emission (black) for 
r ^ 0.15i?vir is about 100 times below that measured in real radio 
halos (red, green). 

By allowing Xp to increase with radius as in |3.3.2| results in an 
increase of the simulated emission for larger radii (figure]?] middle 
panel), but still expectations account for < 10% of the observed 
emission in real radio halos at r ^ 0.15-R v ir- Additional up-scaling 
of the magnetic field (model 3) further increases the size of the 
simulated halos, especially at very large distances from the cluster 
center, where our MHD simulations would predict a steeper scaling 
of the magnetic field with gas density. However this up-scaling only 
mildly increases the level of the expected emission at intermediate 
distances, r w 0.1 — 0.3i? v ir, where the radio brightness of radio 
halos is constrained by present observations. 

In a recent paper Pfrommer et al. ( 2008 1 claim that, accord- 
ing to cosmological simulations, the combined synchrotron emis- 
sion from secondary electrons and primary electrons accelerated 
at large scale shocks may produce diffuse emission with a fairly 
broad spatial distribution. Based on our results, the contribution 
from primary electrons (even assuming Pfrommer et al. (2008 1 re- 
sults) is not expected to solve the discrepancy between models and 
observations in Fig. [7] Indeed, our model 3 is thought to mimic 
the secondary-generated emission in Pfro mmer et al.| (2008| > and, 
based on their Figure 9, the contributions from primary electrons 
at r w 0.2i? v ir is (at best) comparable with that of the secondary 
electrons leaving expectations well below observation^] 

In principle, for each cluster, it would be possible to allow the 
energy density of CRp to further increase at large distances from 
the cluster center and find an ad hoc spatial distribution of CRp that 
allows for matching the brightness profiles of radio halos. However 
this would imply the untenable scenario in which CRp store very 
large energy budget outside the cluster core, for example in the case 
of the Coma halo the energy density of CRp at r fa 0.2 — 0.37? v ir 
should be comparable to that of the thermal ICM (Donnert et al. 
|20T0] >. 

4.3 Morphology of the radio emission 

A complementary approach to compare expectations from simu- 
lations with observations is to derive a point-to-point plot between 
radio and X-ray brightness. A number of clusters indeed show mor- 
phological correlations between X-ray and radio surface brightness 
( |Govon i et al. 2001). This could be preferable to a comparison of 
the radial profiles because morphological distortions and imperfect 
alignment of the emission centers may introduce errors in the ob- 
served profiles ( Govoni et al. 200~T). 

In Figure [8] we show the luminosity in 166 kpc sized patches of 
our synthetic radio maps of the four biggest simulated clusters in 

2 it is more difficult to evaluate the ratio of secondary to primary electrons 
at r ~ 0.3-Rvir from Figure 9 in Pfrommer et al. due to the presence of a 
very bright shock-like spot southest of the center of their simulated cluster 
(see their Figures 7 and 8) that affects the azimuthal brightness profile, but 
that is not indicative of an upturn in the distribution of the truly diffuse 
halo-emission 
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Figure 7. Normalized radial profiles of radio emission from 17 simulated clusters (black). In the left panel the CRe population was modelled using a constant 
Xp. The middle one uses the radius dependent CR to thermal energy fraction (fig. [4) adopted from Pfrommer et al. 1 2008). Additionally we modified (fig. [5) 
the magnetic field in the right panel to be <x ^/p (right panel). We ad observations of A2744, A2319, A545, A2163 and A2255 from (Cassano et al.|2007| as 
red curves with error bars, Coma (green) is taken from Deiss et al. ( 1997) . 
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Figure 8. Radio vers. X-ray surface brightness in patches with 166 kpc side length, with the radio flux normalised to Coma. We also plot a fit to observations 
from |Govoni et al. 1(200 j) (red, b = 0.64 for COMA and A2255.A2319 dotted). The grey box indicates the observed region from |Govoni et al.H200l) . 



Cluster 


Ml 


M2 


M3 





1.09 


1.04 


0.87 


PERSEUS 


1.33 


1.10 


0.90 


COMA 


1.54 


1.18 


0.98 


4 


1.21 


1.01 


0.94 



Table 2. Slopes of the Radio vers. X-ray surface brightness correlation from 
patches for our three models and the four largest clusters. Govoni et al.| 
(20T7l) ; |Rretti et al.|p00T) find slopes of Coma & A2163: b=0.64, A2391: 
b=0.98 and, A2255: b=0.82 



our sample as a function of the same patches in the X-rays. The 
three panels again show constant (model 1, left), increased cosmic 
ray scaling (model 2, middle) and cosmic rays scaled combined 
with magnetic fields scaled (model 3, right). We include the ob- 
served scaling for three clusters (A2255, Coma, A23 19) from |Gov^| 
oni et al. ( 2001 ) (dashed lines) and add the correlations found from 
best fits to the simulated data (dashed dotted lines). In table[2]we in- 
clude the slopes from a best fit model to the simulations inside the 
observed region. A comparison with the values found by |Govoni| 
|eTaLl ( |200l) ; |Feretti et al-lpOOT) , (Coma, A2163: b=0.64, A2391: 
b=0.98, A2255: b=0.82) yields that only the third model may fit 
two of the three observations. Still no model approaches the flat 
value found for the Coma cluster. 

We find that none of the models fit the slopes of any of the observed 
scalings over the whole range, with the radio brightness decreas- 
ing too quickly with respect to the X-ray brightness. This reflects 
the point raised in the previous Section, i.e. that the slope of the 



radial distribution of the synchrotron emission in our secondary- 
generated radio halos is too steep. 

In particular, the adoption of the cosmic ray scaling function 
(models 2 and 3) causes the emission to decrease in the innermost 
five patches, which can be seen in the plots as a bending at high 
brightness. The magnetic field scaling (model 3) adds more power 
to the outmost regions, resulting in a flattening of the correlation, 
but still the obtained slope is steeper than found in observations. 



4.4 Scaling Relations 

There are several observed correlations for radio halos that relate 
thermal and non-thermal properties of the IGM: those between the 
radio power at 1.4 GHz, Pi a, and the X-ray luminosity, Lx, tem- 
perature, and cluster mass (Liang et al.|2000| [Govoni et al.|200i] 
|Cassano et a l. 2006). In addition, by making use of a sample of 
14 giant radio halos, Cassano et al. (2007) found new scaling rela- 
tions that connect the radio power, Pr, of halos to the size of the 
emitting region, Rh (see also Murgia et al. 2009 1, and to the total 
cluster mass within Rh, Mh; a geometrical scaling was also found 
between Mh and Rh- The observed scalings from Cassan o et al.| 
((2007) are : 



Pr 
Pb 



T\/r D 2.17±0.19 

Mh oc H n 



oc Hjj 
oc Ml 



(19) 
(20) 
(21) 



Specifically, Mh was computed from X-ray observations under 
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the assumption of hydrostatic equilibrium and spherical symmetry. 
This procedure may lead to errors as large as 40% in mass ( |Ra- 
sia et al. 2006) which are expected to be not dependent on cluster 
mass, so that these errors might introduce considerable scatter with- 
out affecting the real trend of the correlation. Rn was measured on 
the radio images, Ru = V-Rmin X Pmax, where R m in and P ma x 
are the minimum and maximum radii measured on the 3<r radio 
isophotes. We stress that Rn provides a simple, but viable estimate 
of the physical size of radio halos, indeed a one-to-one correlation 
has been found between Rh and the size containing the 85% of the 
radio halo flux, Rss, derived from the observed brightness profiles 
of halos ( Cassano et al. 2007}. 

A scaling was also found between the size of radio halos and 
the virial radius of clusters, Ph oc R^'® 3±0 5 ° ^Cassano et al 



2007). Given that massive clusters are almost self similar (e.g. 



Rosati et al. 2002) one might have expected that Ph scales with 



R vl r and that the radial profiles of the radio emission are self- 
similar. On the contrary, this result proves that self-similarity is 
broken in the case of the non-thermal cluster components, as first 
noted by (Kempner & Sarazi n~]2001[ >. 

As the synchrotron power depends on both magnetic field scaling 
and CRe scaling with density, it is unclear what is responsible for 
the break in the observed properties. On the other hand we know 
from previous work (Donnert et al. 2008) that the magnetic field 
scaling (with temperature or mass) flattens out for the largest clus- 
ters in our simulation, that would imply an expected break of self 
similarity in the thermal vs non-thermal properties of our simulated 
clusters. 

|Cassano eta l. ( 2007 1 showed that all the correlations explored 
so far for radio halos can be derived by combining the Ph — P„ and 
Pi. 4 — Rh scalings. This suggests that there are two main scaling 
relations that carry out the leading information on the physics of 
the non-thermal components in galaxy clusters. 

In what follows we shall investigate whether the properties 
of our simulated secondary-radio halos are consistent with the ob- 
served scalings. 



4.4.1 Mass vs. Size 

As a first step, before comparing the observed scalings with those 
derived for our simulated clusters, we check whether our clusters 
inherit the same mass distribution of real clusters with radio halos. 
We compare the observed and simulated scaling between the total 
mass inside Rh (Mss) and R$ 5 , that provides a geometrical scaling 
on the halo-region. Therefore we plot in Figure [5] Mn versus Pgs 
for simulated (open triangles) and observed (diamonds) clusters, 
together with the best fit power law to the observed scaling from 
|Cassano et aTlpOOT) . 

The three models define different values of Pgs and conse- 
quently different volumes where the scaling can be tested. We find 
that in all cases simulated clusters lie on the thermal scaling de- 
scribed by observed clusters, although, as expected (Sects [4~2]|4-.3^ 
simulations populate a region in the Mn-Rss diagram with smaller 
values of Rss- 



4.4.2 The Size vs. Size relation 

As stated in Section |4~4"1 observations of clusters with radio halos 
show that the size of halos scales not linearly with the cluster virial 
radius, suggesting that the non-thermal component in clusters is not 
self-similar. 



In figure [10] we plot the radius containing 85% of the clusters 
emission over virial radius for all three models. We include data 
from |Cassano et aL[ j2"007) and the fits to the simulated clusters- 
distribution obtained for models 1-3 at z = 0. We find a correlation 
between Pss and P v j r for our simulated hadronic-halos. Results 
suggest that self-similarity is preserved in the non-thermal compo- 
nents, as the increase of the halos's radius is roughly proportional 
to the virial radius of the hosting clusters; the slopes of the correla- 
tions are 6f = 0.8, 1.2 and 1.3 for models 1, 2 and 3, respectively. 
This is not in line with observations: the expected correlations are 
flatter than the observed one and we predict halos systematically 
smaller than the observed ones. In contrast this is expected con- 
sidering results reported in section [4~2| and confirms that it is chal- 
lenging to reproduce the extension of the observed radio halos with 
hadronic models, even by adopting a profile of the magnetic field 
that is flatter than that from our MHD simulations and by assuming 
a flat profile of the spatial distribution of CRp (model 3). 
|Donnert eTa l. (2010) have shown that matching the radio emission 
of the Coma halo at distance ~ 0.2 — 0.3 R V i r with hadronic mod- 
els (by further inceasing the CRp energy content at larger radii, 
see also section [4~2] l would require the energy content of CRp to be 
roughly similar to the content of the thermal ICM at these distances 
(at least when constraints on the magnetic field from RM observa- 
tions of the Coma cluster are used). The scalings found in figure [TO] 
make the situation possibly more challenging, because they show 
that even more energy in the form of CRp would be required in the 
case of more massive clusters where the differences, between Pgs 
of our hadronic-halos and that of the observed ones, are larger. 

The observed clusters are in a range of redshifts, z — — 
—0.4, while hadronic-halos are extracted from the simulations at 
z — 0. For this reason in figure [Tofc we also plot the distribution 
of simulated clusters obtained at z = 0.4 in the case of model 
1 (marked green). We find that the correlation between Rss-R v ir 
does not show significant change in slope at higher redshift. 

4.4.3 The X-ray Luminosity vs. Radio power relation and the 
evolution of radio halos 

Radio halos follow a correlation between the monochromatic radio 
luminosity at 1.4 GHz, Pi. 4, and the X-ray luminosity of the host- 



ing clusters, Lx (e.g. Liang et al. 


2000 


EnBlin & Rottgering 2002; 


Bacchi et al. 2003 , Cassano et al. 


2006 


1. Recent radio observations 



of a statistical sample of X-ray selected galaxy clusters, the "GMRT 
radio halo survey", ( |Venturi et al.|2007| [2008] ) allow to study the 
distribution of clusters in the Pi. 4 — Lx diagram. These obser- 
vations suggest that the distribution of clusters in the Pi. 4 — Lx 
diagram is bi-modal: radio-halo clusters trace the Pi. 4 — Lx cor- 
relation, while the majority of clusters are found "radio quiet" with 
the limits to their radio luminosities about 10 times smaller than the 
radio luminosities of halos. 

In order to investigate the behaviour of our simulated clusters in 
the Pi. 4 — Lx diagram, in Fig[TT|we plot our simulated clusters 
together with observed clusters (from Ve nturi et al.| ([2008 1). The 
synthetic radio luminosities of our simulated clusters are scaled in 
order to have the simulated Coma cluster matching the observed 
one. 

According to secondary models a correlation between radio 
luminosity and cluster X-ray luminosity (or temperature) is ex- 
pected (e.g.|Dolag & EnMin|2000[|Miniati|200r)|Dolag et al.|2005[ 



|Pfrommer et al.|2008| l. We qualitatively confirm these expectations 
and in all three models find that the largest simulated clusters would 
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Figure 9. Total gravitational mass inside the radio emitting region over radius of the same region. We plot the correlation for all three models, constant and 
varying CR fraction and varying CR fraction and upscaled magnetic field (from left to right). We also plot observations of 14 galaxy clusters from iCassano 
|et al.|2 007 1 and the best fit correlations. Due to uncertainties in the mass estimation from X-ray data the observations may show systematic errors. For both, 
simulations and observations the coma cluster is marked in blue. 
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Figure 10. Radius of the radio emitting region as function of the virial radius of the cluster. We plot, left t o right, the correla tion for all three models and 
include the best fit (x 2 = 0.54, 0.29, 0.86, respectively) . We also plot observations of 14 galaxy clusters from {Cassano et al.|20f?7| and the best fit correlation. 
For both, simulations and observations the Coma cluster is marked in blue. For model 1 we include the include the correlation at redshift 0.4. 



naturally approach the observed correlation. In all three models the 
smallest systems are significantly more scattered in radio power 
than the largest clusters. This is expected as the magnetic field in 
the central regions of our simulated clusters (where most of the 
synchrotron emission is generated) is found to be tighly correlated 
with cluster thermal properties only in the case of massive clusters, 
while only a steep trend is found for smaller systems (figure 9 in 
|Donnert et al.|2008| >. 

The most relevant difference with respect to observations is 
that according to all three models (as for every secondary model) 
the synchrotron luminosity of the massive simulated clusters is 
equivalent to that of typical radio halos, at least if the radio lu- 
minosity of the simulated Coma cluster is normalised to that of 
the real Coma halo. This is inconsistent with observations which, 
on the other hand, found radio halos in only about 1 /3 of massive 
clusters. 

Most important, no radio bi-modality is expected in our simu- 
lated secondary-halos. Hadronic halos in simulated massive sys- 
tems would follow a tight correlation, while those in less massive 
systems would be more broadly distributed. 

One possibility to reconcile the hadronic scenario with the observed 
halo-merger connection and with the bi-modal distribution of clus- 
ters in the Pi. 4 — Lx diagram is to admit that the observed bi- 
modality is driven by the amplification and dissipation of the mag- 
netic field in the merging and post-merging phase, respectively 
( |Brunetti & Lazarian|[2007l |Pfrommer et al.|[2u08l |Kushnir et al.| 
|2009| >. However, Burnetti et al.(2009) have shown that the degree of 
amplification/dissipation of the magnetic field and the time-scale 
of this process that would be necessary to explain observations are 



difficult to reconcile with the observed properties of magnetic fields 
in the ICM (namely the field intensity and coherence scales from 
Rotation Measurements), and appear also disfavoured by energetic 
arguments. 

The amplification of the magnetic field during cluster mergers is 
followed by our MHD cosmological simulations that offer a com- 
plementary approach to highlight this issue. To show the effect of 
magnetic field evolution and investigate the bimodality we plot the 
time evolution of the simulated radio haloes in figure^] Shown is 
the radio luminosity over X-ray luminosity of the simulated sample 
for redshift z < 0.48 (big triangles: z=0, small triangles and black 
line: earlier redshifts); also we show the radio power vers, temper- 
ature correlation and evolution in the Appendix. No hint of a bi- 
modality is found, simply because magnetic field amplification in 
massive systems is a gradual process that happens in a time-scale 
comparable with the life-time of clusters themselves. The evolution 
of magnetic field is reflected in the broad/scattered distribution of 
clusters in the Pi, 4 — Lx diagram, especially in the case of smaller 
systems. These smallest clusters approach the lower end of the cor- 
relation with rather large scatter because the magnetic field is not 
saturated yet in their central regions and even smaller mergers yield 
a significant field amplification resulting in an increased radio lu- 
minosity. Although one would expect this behavior to be slightly 
dependent on numerical resolution, increasing the numerical res- 
olution would make the situation even more stringent. Resolving 
smaller gas motions leads to increase the amplification of the mag- 
netic field, specially at early times. Having a larger magnetic field 
the changes in the magnetic field due to merger activity will be sup- 
pressed as the magnetic field is already closer to saturation effects. 
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Figure 11. Radio power per frequency at 1.4 GHz over X-ray luminosity in erg/s from our simulated clusters using all three models (triangles, ltr.: constant 
fraction, scaled fraction, scaled fraction and scaled field). Observed scalings by Cassano et al. 1 2007 1 (diamonds, z > 0.2: red diamonds) and non-detections 
in red I Venturi et al. 2007 2008). For both, simulations and observations the Coma cluster is marked blue. 
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Figure 12. Radio power per frequency at 1.4 GHz over X-ray luminosity 
in erg/s from our simulated clusters using the constant model. Observed 
scalings by Cassano et al. (2007) (diamonds, z > 0.2: red diamonds) and 
non-detections in red (Venturi et al.|2 007 2008 1. We also include the time 
evolution of clusters for z < 0.48 (black lines) 

Therefore we would expect the clusters (especially the smaller, less 
resolved systems) to evolve even more along the correlation, and 
the spread around it would be further decreased. 



5 7-RAY EMISSION FROM SIMULATED CLUSTERS 

We use the formalism in section |3.2| to compute predictions for 
the 7-ray luminosity of the simulated clusters according to all 
three models. Of interest are the two energy bands of Cherenkov 
telescopes, E > 0.1 TeV, and FERMI/EGRET telescopes, 
£>0.1GeV. 

In table [3] j4]l we present fluxes for the simulated sample in the 
VERITAS (FERMI) energy band. None of the clusters is in the 
observable range of the VERITAS experiment, however the largest 
ones (0, Coma, Virgo, Perseus, Centaurus) have a chance to be 
detected by FERMI in next years, at least for models with radially 
increasing CRp normalisation (model 2 and 3). 

In the left panel of figures [T3] and [14] we plot the simulated 
clusters in 7-ray vers, radio luminosity at 1.4 GHz for the 3 
models. Following the previous Section we normalise the radio 
luminosity of the sample so the simulated Coma cluster fits the 



Figure 13. 7-ray luminosity (E > 0.1 Tev, e.g. VERITAS) over radio lu- 
minosity (at 1 .4 GHz, left) and bolometric X-Ray luminosity of all cluster 
at redshift z=0.We plotted the three models as different symbols (diamonds 
- model 1, triangles - model 2, boxes - model 3) and include upper limits 
from|Perkins 12008 1; Per kins & the VERITAS Collaborauon1(2056) in blue. 
The Coma cluster is marked green. 



observed emission from the Coma halo. We also report available 
limits from VERITAS observations jPerkins|2008l |Perkins & The] 
I VERITAS Collaboratioril[2006l > (figure [13). Further we include 
(figur e|l4| n ewest preliminary results from the FERMI experiment 
(Mori 2009) on the upper limit of the 7-ray flux from the Coma 
cluster. A scaling between 7-ray and radio luminosity is expected 
for the 3 models because p-p collisions generate both secondary 
electrons and neutral pions, the scaling is sub-linear since the 
radio luminosity is further boosted by the increase of the magnetic 
field in more massive clusters. 

In the right panel of figure [T3| and 1 1 4| we plot the emission 
of our clusters in 7-ray vers. X-ray luminosity for the 3 models. 
Also in this case we show two 7-ray upper limits from VERITAS 
( |Poulain et al.|1992^ and preliminary FERMI jMori|2009| obser- 
vations respectively. A quasi-linear correlation is predicted, and is 
found less scattered than that between 7-ray and radio luminosities 
since it compares two purely thermal quantities (also in the case of 
models 2 and 3 CRp are scaled with thermal energy according with 
the profile in Figure[4](. 

The FERMI results are still consistent with the models for CRp pre- 
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Figure 14. 7-ray luminosity (E > 0.1 Gev, e.g. FERMI, EGRET) over 
radio luminosity (at 1 .4 GHz, left) and bolometric X-Ray luminosity of all 
cluster at redshift z=0.We plotted the three models as different symbols (di- 
amonds - model 1 , triangles - model 2, boxes - model 3) . We add prelimi- 
nary results from the first year FERMI data I Mori 2009 1 in red. The Coma 
cluster is marked green. 
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Table 3. Fluxes (7/cm 2 /s) in the VERITAS energy range (E > 100 GeV) 
from 1 ilvir for all Clusters from the simulations using the three different 
models. 



sented in this paper. As shown by Donnert et al. (2010) the expected 
7-ray flux increases substantially if we consider a spatial CR pro- 
ton distribution so the radio emission actually fits the observations 
in morphology. We therefore expect a rejection or confirmation of 
our models in the very near future. 



6 CONCLUSIONS 

We use a constrained, cosmological MHD SPH simulation with 
a semianalytic model for galactic magnetic outflows, to obtain a 
sample of 16 galaxy clusters with thermal properties similar to 
clusters in the Local Universe. Further we assume 3 different mod- 
els for secondary cosmic rays motivated by simulations (Pfrommer 
|et al.| [2007l using the proper high energy approximation for the 
pion cross-section. In the first model we assume a constant CRp 
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Table 4. Fluxes (7/cm 2 /s) in the EGRET/FERMI energy range (E > 
0.1 TeV) from 1 R v i r for all Clusters from the simulations using the three 
different models. 



normalisation relative to the thermal density and the simulated 
magnetic field. In the second model we keep the magnetic field 
and introduce a radius dependent CRp normalisation infered from 
non-radiative simulations by |Pfrommer et al.| |2007 ). The third 
model uses the same CRs and flatten the simulated magnetic field 
to be B oc ^fp. 

Although our simulations do not include a (internally) self- 
consistent treatment of CRp as done in other simulations (e.g. 
|Pfrommer et al.||2007) , contrary to previous work they allow to 
properly simulate the properties of the magnetic field in the ICM 
which is important for modeling the cluster-scale synchrotron emis- 
sion. 

For the first time, we carry out a detailed comparison between 
the observed properties of giant radio haloes and those of simulated 
halos according to secondary models and under different assump- 
tions for the spatial distribution of CRp. 

In an earlier paper we presented a detailed comparison of 
the simulated Coma cluster (Donnert et al. 2010). In this work 
we focus on global sample properties and compare with recent 
observations. 

In particular, as a first step, we show that : 

• The radial profiles of Faraday rotation of the median of our 
cluster sample is in line with that obtained from a number of ob- 
servations of different clusters. This confirms that the properties of 
the clusters's magnetic fields in our simulations are rather similar 
to the observed ones. 

• The normalised radial profiles of the radio emission at 1.4 
GHz of our simulated hadronic-halos show a deficit at radii ^ 
0.1r v ir with respect to the synchrotron profiles observed for a sam- 
ple of well studied radio halos. This is in line with previous claims 
based on semi-analytic calculations in the context of the hadronic 
model (Brunetti et al. 2004). In addition our results show that, even 
by assuming a flat profile for both the magnetic field and CRp spa- 
tial distributions (Model 3), secondary electrons may account for 
less than about 10-15 % at radii 0.15 - 0.3 r viv . 

• A Point to Point comparison of radio vs. X-ray emission, ob- 
tained for the 4 largest clusters of our sample, confirms the results 
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obtained for the profiles showing that the radial distribution of radio 
emission is too steep. All three models do not fit the observations 
over the whole range. Furthermore an excess in radio emission of 
the innermost patches suggests that haloes from the simulation are 
to centrally peaked. 

As a second step we compare scaling relations obtained for the 
hadronic-halos in our simulated cluster sample with those given in 
( |Cassano et al.|20 07 1 that are obtained from a sample of observed 
radio halos. We find that : 

• The geometrical correlation between the radius of radio halos 
and the cluster-mass contained within this radius is well reproduced 
by all three models. Due to the expected self-similarity of cluster 
in thermal properties this result implies that -at least- the simulated 
and observed clusters share similar physical properties. 

• A quasi self-similar behaviour is found for the non-thermal 
properties of our simulated clusters. In particular the radius of our 
hadronic halos is found to scale (approximately) with the virial ra- 
dius of the simulated clusters. This is contrary to observations that 
found a steeper correlation between halo-radius and virial radius of 
the hosting clusters and implies that our simulated halos are sys- 
tematically smaller than the observed ones. 

• A correlation between the monochromatic luminosity of our 
hadronic halos and the X-ray luminosity of the simulated hosting- 
clusters is found. As soon as the population of our hadronic halos is 
normalised, by scaling the radio luminosity of the simulated Coma 
halo with that of the observed one, the correlation is similar to that 
observed for radio halos. However, since at this point all the sim- 
ulated clusters show radio emission at the level of the observed 
halos, we find that the cluster radio bi-modality, observed for X- 
ray selected clusters, cannot be reproduced. A radio bi-modality 
in the radio - X-ray diagram would require a fast (< Gyr) evo- 
lution of the radio luminosity in connection with cluster mergers, 
on the other hand, we find that the time evolution of our simulated 
massive-clusters in this diagram happens on cosmological, long, 
time-scales. Finally, we show that once the radio - X-ray correla- 
tion is approached at low X-ray luminosities, clusters follow the 
correlation closely due to the saturation of the magnetic field. 

As a final point we calculate the 7-ray emission from our sim- 
ulated clusters once the radio luminosity of the simulated Coma 
halo is anchored to that of the observed one (essentially by scaling 
the number density of CRp in simulated clusters). We find that : 

• The 7-ray emission expected from our simulated clusters is 
well below the sensitivity of present Cherenkov Arrays, e.g. the 
VERITAS experiment, for all the adopted models for hadronic ha- 
los. The 7-ray fluxes at > 100 MeV expected from our simulated 
clusters would allow for a marginal detection by the FERMI tele- 
scope in next years, at least by assuming Models 2 and 3. 

• The integrated 7-ray flux from our simulated clusters is ex- 
pected to scale with their radio emission, although with rather large 
scatter. On the other hand, a tight correlation is found between 7- 
ray and X-ray fluxes, in which case (due to the scalings between 
thermal and non-thermal CRp in Models 1-3) the correlation is es- 
sentially driven by the density of the thermal gas in our simulated 
clusters. 

Considering all these results, as well as the outcome from our 
earlier work on the Coma cluster (Donnert et al. 2010), we con- 
clude that hadronic models alone are not able to explain the ob- 
served properties of giant radio haloes, in terms of their radial ex- 
tension, observed clusters's radio bi-modality, scaling relations and 



spectral properties. Therefore we conclude that radio emission ob- 
served in galaxy clusters in form of giant radio halos are a powerful 
tool to infer the amount of cosmic ray protons within galaxy clus- 
ter, confirming previous results that the energy content of cosmic 
ray protons in clusters can not exceed percent level (e.g. |Brunetti| 
|et al.|2007||Churazov et al.|2008) . 

The problem of halo's extension could be alleviated in ex- 
tended hadronic models, where the contribution from primary 
(shock accelerated) electrons at the cluster outskirts is combined 
in the simulations with that from secondary electrons ( |Pfrommerj 
et al. (2008)). However, as shown in Sect. |4.2| the problem of the 
halo profiles arises at distances 0.1 - 0.3 R V1I from cluster centers 
where, based on the same simulations, the contribution to the dif- 
fuse synchrotron emission from these shock accelerated electrons 
is not yet dominant. Further we would like to note that the detailed 
morphology of the radio emission caused by electrons injected at 
shocks (e.g. |Roettiger et al,|1999 10 will be significantly different 
from the morphology of giant radio halos and therefore is in general 
unlikely able to mimic large scale, cluster centric radio emission. 
Thus additional physical processes are required to explain observa- 
tions. 

The problem of the radio bi-modality in massive clusters 
can be alleviated by assuming that MHD turbulence (and the rms 
magnetic field) decays as soon as clusters approach a relaxed state 
after a major merger ( |Brunetti et al.|2009|[Kushnir et al.|2009) . In 
this case it might also be thought that CRp in these relaxed clusters 
would undergo less scattering on magnetic field irregularities 
escaping from the cluster volume and reducing the source of 
secondary electrons. All these processes are not properly included 
in our MHD cosmological simulations. However, Brune ttTet al.| 
(2009) have shown that the fast decay of the rms field in galaxy 
clusters, that is necessary to explain the observed bi-modality, 
would imply a rather unphysical situation where the magnetic 
field power spectrum peaks at the smaller scales, and leads to the 
consequence that an extremely large flux of energy in clusters goes 
into magnetic field amplification/dissipation. In addition, because 
turbulent cascade should start to decay on largest scales, which are 
resolved by Faraday rotation measurements, this scenario would 
predict a bi-modality in rotation measures and depolarisation in 
largest clusters, which is not observed (Govoni et al. 2010). 

The discovery of giant radio halos with very steep spectrum 
(a ~ 1.8 — 2) in merging clusters proves that inefficient par- 
ticle acceleration mechanisms are responsible for the origin of 
these sources and, based on simple energy arguments, disfavours 
hadronic models ( |Brunetti et al.||2008| >. Observations of the best 
studied halo, in the Coma cluster, have shown that its spectrum 
steepens at higher frequencies (Thier bach et al.|20 03 ), due to the 
competition between energy losses and acceleration of the emit- 
ting particles, in which case the particle acceleration time-scale 
would be of the order of 0.1 Gyr. More recently we have shown 
that the steepening cannot be a result of the inverse Compton (SZ) 
decrement (Donnert et al. 2010), confirming these previous finding 
and supporting a scenario of inefficient particle acceleration mech- 
anisms at the origin of the halo. 

Particle acceleration due to micro-turbulence in merging galaxy 
clusters has been proposed for the origin of radio halos (e.g. 
Brunetti et al. 2001; Petrosian 2001). In this case gentle particle 



3 These are base on Eularian simulation, which are most appropriate to 
describe the morphology of shocks outside the cluster core. 
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acceleration mechanisms would generate radio halos in connection 
with (massive) cluster mergers, while the radio emission would de- 
cay as soon as clusters approach a relaxed state, due to dissipation 
of (at least a fraction of) this turbulence and the fast electron radia- 
tive cooling. 

It will therefore be necessary to consider these processes in 
future simulations used to study radio halos. That includes an esti- 
mation of the locally merger injected turbulence as well as a more 
detailed description of CR electron spectra. 
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Figure Al. Bolometric X-Ray luminosity over temperature of all cluster at 
redshift z=0 (diamonds) and z=0.48 (triangles). We also plot the observed 
correlation from |Pratt et al.H2008) . 
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Figure A2. Radio Power at 1.4 Ghz over mass weighted temperature inside 
0.1r v ; r for all clusters. We include various observations from the literature. 
The evolution of the simulated clusters is plotted as lines in black (gray) for 
z < 0.48 (z > 0.48), while the points correspond to z = 0. 



Here we provide additional scaling relations to substantiate the 
scalings shown before. 



Al X-Ray luminosity - Temperature relation 

In figure |A1| we show bolometric X-ray luminosity over mass 
weighted temperature for the simulated cluster sample at redshift 
zero (diamonds) and 0.48 (triangles). We include band-corrected 
the observed relation from Pratt et al. ( 2008 1 (broken line) and a the 
same correlation normalised to the simulated Coma cluster (black 
line). It is well known that non radiative simulations tend to over 
predict the x-Ray luminosity by a fair amount (see for example 
(Borgani et al. 2006 1). In addition the effect of K-correction can be 
seen in comparison with the cluster at high redshift. 
Still, after re-normalizing the x-Ray luminosity there is a reason- 
ably good agreement of with observations We therefore conclude, 
that our simulated clusters sample shows similar thermal properties 
than the observations. 



A2 Radio luminosity - Temperature relation 

In figure |A2| we present the scaling of radio power at 1.4 GHz over 
cluster temperature as obtained from e.g. X-ray observations for 
constant CRp scaling (model 1). We include the evolution of the 
emission with time as lines in black (gray) for z < 0.48 (z > 0.48). 
Further we plot a number of recent observations. The sample fol- 
lows the correlation closely at redshift zero, and only two clus- 
ters show significant deviation from the correlation at temperature 
larger than 5 keV, and only at high redshift. 

In simulations, temperature is a less merger sensitive mass estima- 
tor compared to the X-ray luminosity, in terms of cluster mergers. 
We therefore conclude, that the bimodality observed in large clus- 
ters is not a result of biased mass estimation. 



